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ABSTRACT 

O 1 A thorough analysis of multicolour CCD observations of two modulated RRab-type 

| variables, XY And and UZ Vir is presented. These Blazhko stars show relatively simple 

light-curve modulation with the usual multiplet structures in their Fourier spectra. 
One additional, independent frequency with linear-combination terms of the pulsation 
frequency is also detected in the residual spectrum of each of the two stars. The 
amplitude and phase relations of the triplet components are studied in detail. Most of 
1 the epoch-independent phase differences show a slight, systematic colour dependence, 

0^ . however, these trends have the opposite sign in the two stars. The mean values of 

C*") ■ the global physical parameters and their changes with Blazhko phase are determined 

utilizing the Inverse Photometric Method (IPM). The modulation properties and the 
IPM results are compared for the two variables. The pulsation period of XY And is 
OO ' the shortest when its pulsation amplitude is the highest, while UZ Vir has the longest 

\ pulsation period at this phase of the modulation. Despite this opposite behaviour, the 

■ phase relations of their mean-physical-parameter variations are similar. These results 

are not in accord with the predictions of the Blazhko model of Stothers (2006, ApJ, 
652, 643). 

■ Key words: stars: horizontal branch - stars: variables: other - stars: individual: XY 
5_i " And, UZ Vir - stars: oscillations (including pulsations) - methods: data analysis - 

techniques: photometric 



1 INTRODUCTION 

In spite of the recently increased research activity in the 
field, we still do not have a satisfactory model for the phe- 
nomenon of the light-curve modulation of RR Lyrae stars, 
the so-called Blazhko effect, that would be consistent with 
the wide variety of the ever increasing collection of obser- 
vations. To find such a model, not only the phenomenol- 
ogy of the light-curve variation have to be known; the 
modulational changes in the atmospheric physical param- 
eters of these stars also impose constraints on the poten- 
tial Blazhko models. The brightness variations of Blazhko 



* Based on observations collected with the automatic 60-cm tele- 
scope of Konkoly Observatory, Budapest, Svabhegy 
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RR Lyrae stars can be studied with unprecedented de- 
tail and precision through the data of CoRoT and Kepler 
space telescopes in a single photometric band, and these 
studies led to new explanations of the phenomenon with 
the detection of half-integer frequenci es (jBuchler fc Kollathl 
120121 : iKollath. Molnar fc SzabcTbOllI '). However, to derive 
physical parameters, standard multicolour or spectroscopic 
observations are needed. Since the modulation periods of 
Blazhko stars span a wide range from about one week to 
many hundred days, detailed spectroscopic studies with suf- 
ficiently large-aperture telescopes are rather difficult to at- 
tain. Therefore, moderate-size earth-based multicolour pho- 
tometric telescopes still have an important role in studying 
the Blazhko effect. 

In the context of the Konkoly Bla zhko Survey I and 
II (|jurcsik et al.ll2009l ; ISodor et alj|2012h , we have obtained 
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Table 1. Log of the CCD observations of XY And and UZ Vir 
obtained with the 60-cm automatic telescope of Konkoly Obser- 
vatory. 



Object 


From 
(JD- 


To 
2 450 000) 


Filter 


Nights 


Data 
points 


XY And 


4390 


4857 


V 


64 


2611 








la 


64 


2610 


UZ Vir 


4504 


4978 


B 


72 


1722 








V 


70 


1695 








Ic 


70 


1676 



Table 2. Identification of the variables and their comparison 
stars. 



Object 



2MASS ID 



Comp. star 2MASS ID 



XY And 01264243+3404068 
UZ Vir 13084432+1324084 



01270016+3404307 
13082756+1322403 



multicolour light curves of bright, northern, fundamental- 
mode Blazhko RR Lyrae stars. These data have already 
been successfully utilized to derive the variations in the at- 
mospheric paramet ers of six modulated stars (see sect. 6 in 

ljurcsik et al.|[20l3 and references ther ein). 

According to the explanation of IStothersi (|2006l ). the 
Blazhko effect is 'a direct consequence of a gradual strength- 
ening and weakening of turbulent convection in the stel- 
lar envelope'. Though this m odel has been s eriously criti- 
cized on the oretical grounds by Kovacsl (l20Qgh ; ISmolec et alj 
( |201lf ) and iMolnar. Kollath fc Szabd l|2012l ). it is also im- 
portant to confront the predictions of the Stothers-model 
with observational results. The model predicts that the 
phase relation between the amplitude and period changes 
depends on the mean physical properties of the star and 
it is also in direct connection with the phase relations of 
the physical parameter variations during the Blazhko cy- 
cle. The present paper aims to test if there is any connec- 
tion between the phase relation of the amplitude and pe- 
riod (phase) variations and the mean physical parameters 
an d their variations dur ing the Blazhko cycle as proposed 
by IStothersi |2006j . l201ll ). Therefore, we have selected two 
stars to study their Blazhko modulation in detail, XY And 
(02000 = 01 h 26 m 42?43, 5 2 ooo = +34°04 , 06'.'9, P = 0.3987 d, 
Pmod = 41.4 d) and UZ Vir (a 20 oo = 13 h 08 m 44=32, 5 2000 = 
+ 13°24'08'.'4, P = 0.4593 d, P mod = 68.2 d), which show dif- 
ferent phase relations between their amplitude and period 
(phase) variations. XY And and UZ Vir were extensively 
observed in the cou rse of the Konkoly Blazhko Survey I 
ijjurcsik et al.ll2009h . but the data and the results of the 
light-curve analysis have not been published yet. 



2 OBSERVATIONS 

Multicolour CCD observations of XY And and UZ Vir were 
obtained between February 2008 and May 2009, and be- 
tween October 2007 and January 2009, respectively, with 
the 60-cm automatic telescope of Konkoly Observatory at 
Svabhegy, Budapest, equipped with a 750 x 1100 Wright 
Instruments CCD camera. In the case of UZ Vir, Johnson- 



Cousins BVIc filters were used, while for the ~ 0.7 mag 
fainter XY And, the observations were made only in the V 
and 7c bands in order to assure a reasonably short duty 
cycle. In each band, the total numbers of the observations 
are about 2600 and 1700 for XY And and UZ Vir, respec- 
tively (see Table [1] for the log of observations.) Standard 
CCD calibration and aperture photometry were performed 
on all the object frames using IRAlQ packages. The rela- 
tive magnitudes of the variables were measured against the 
nearby comparison stars identified in Table [2] The relative 
magnitudes were transformed to the Johnson-Cousins sys- 
tem utilizing standard procedures. 

Fig. [TJshows the light curves of the variables folded with 
the pulsation period at selected Blazhko phases. 

The light-curve data of XY And and UZ Vir are avail- 
able as supporting information with the online version of 
this paper. Tables lATrlXBI in Appendix \K\ give information 
regarding the form and content of the data. 



3 FREQUENCY ANALYSIS 

The mathematical descripti on of the light cu rves is sought by 
Discrete Fourier Transfor m (iDeemind 19751 ) as implemented 
in the mufran package (|Kollathl 1990l ), the gnuplo t utili- 
tiefi and a non-linear fitting algorithm (|S6dorll2012i ). 

A peak in the Fourier spectrum is accepted as an inde- 
pendent frequency of the star and is included into the solu- 
tion if its amplitude exceeds 3.5<r in at least two photometric 
bands (a is denned as the local mean level of the residual 
spectrum). For linear combination terms of the independent 
frequencies, the selection criterion is that its amplitude ex- 
ceeds 2.5<r in at least two bands or 3<r in one band. 



3.1 Pulsation and modulation components 

The top panels of Fig. [2] show the spectra of the V light 
curves of XY And and UZ Vir. The bottom panels display 
the residual spectra after prewhitening for the pulsation fre- 
quency and its harmonics (kfo). The window functions are 
plotted in the inserts. Side-peaks belonging to the modula- 
tion triplet frequencies (kfo ± / m ) with significantly differ- 
ent amplitudes of the two components are clearly identified 
in the residual spectrum of the variables. After prewhiten- 
ing for the pulsation frequency, its harmonics, the modula- 
tion frequency and their identified linear-combination terms, 
no further triplet (multiplet) structure, indicative of a sec- 
ondary modulation, is found around the pulsation frequency 
and its harmonics in any of the stars. 



3.2 Detection of additional frequencies 

Fig. [3] shows the residual spectra of XY And and UZ Vir in 
the observed bands after the removal of the detected multi- 
plets. Both spectra have a regular pattern of peaks with sep- 
arations of the pulsation frequency, /o- The frequency values, 

1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

2 http://www.gnuplot.info/ 
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Figure 1. Folded V light curves of XY And and UZ Vir in selected Blazhko phase bins. The gray dots indicate all the measurements. 
In the top-left corners, the Blazhko-phase ranges of the highlighted data are given, according to the Blazhko cphemerides, Eqs.[2]and 
[4] The rising branch of XY And shows a significant phase variation (~0.1 pulsation phase) while the fix point on the rising branch of 
UZ Vir indicates that amplitude modulation dominates its light-curve variation. The phase of the maximum brightness varies, however, 
within about the same, ±0.02 pulsation-phase range in both stars, as shown in Fig. [4] 



XY And 

1f 2f 3f 4f 5f 




-0.10 -0.05 0.00 0.05 0.10 

Frequency [c/d] 



-0.10 -0.05 0.00 0.05 0.10 

Frequency [c/d] 



Figure 2. Fourier amplitude spectra of the V light curves of XY And and UZ Vir (top panels). The residual spectra are plotted in the 
bottom panels around the first four harmonics of the pulsation frequency after the removal of the detected pulsation components (fc/n, 
k = {1, .... 13} for XY And and k = {1, 25} for UZ Vir). The spectrum segments are artificially shifted for clarity. Peaks belonging to 
the modulation (kfo ± / m ) are clearly identified. The window functions, shown in the inserts, demonstrate that no high-amplitude alias 
frequencies are present besides the daily and yearly aliases. 
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Figure 3. Residual amplitude spectra of XY And and UZ Vir in 
the observed bands after the removal of the pulsation frequency, 
its harmonics, the modulation frequency and their linear combi- 
nations. Dashed lines show 3.5 times the mean level of the spectra 
(3.5<r). A regular pattern is evident in the spectra of both stars: 
peaks are present at 5.21, 7.72, 10.23 and 12.74cd~ 1 , and at 2.89, 
5.06, 7.24 and 9.42 cd" 1 in the case of XY And and UZ Vir, re- 
spectively. These frequency series (indicated by vertical dotted 
lines) are separated by the pulsation frequency. 



the fo/f frequency ratios, the amplitudes and their signifi- 
cance are given in Table [3] Most probably, one of these fre- 
quencies is an additional (radial or non-radial) mode, while 
the other peaks are linear combinations of this mode and the 
pulsation frequency. For UZ Vir, ±lcd _1 alias solutions for 
the identification of the 'base' frequency of the additional 
frequency series cannot be excluded either. 



Table 3. The sequences of identified additional frequencies in 
XY And and UZ Vir. The ratios with the pulsation frequency 
(/o//j')> the amplitudes and the significance are given. For UZ Vir, 
the ±1 cd~ 1 alias frequencies are listed as well. 









XY And 








freq. 


fo/fi 


A(B) 


(7 


A(V) 


a 


Mi) 


IT 


(cd- 1 ) 




(mag) 




(mag) 




(mag) 




/{ 5.2122 


0.481 






0.0028 


3.3 


0.0013 


2.2 


f 2 7.7202 


0.325 






0.0032 


3.6 


0.0019 


3.5 


f' z 10.2282 


0.245 






0.0026 


3.5 


0.0011 


2.2 


f\ 12.7362 


0.197 






0.0022 


3.1 


0.0011 


2.3 


UZ Vir 


f[ 2.8880 


0.754 


0.0033 


3.8 


0.0027 


3.9 


0.0019 


2.4 


f! z 5.0648 


0.430 


0.0023 


3.0 


0.0029 


1.5 


0.0025 


3.9 


/| 7.2416 


0.301 


0.0032 


4.2 


0.0025 


4.2 


0.0021 


3.7 


f' A 9.4184 


0.231 


0.0019 


2.8 


0.0017 


3.0 


0.0018 


3.7 


f[~ 1.8880 


1.153 


0.0032 


3.9 


0.0022 


3.3 


0.0015 


2.1 


f!f 4.0648 


0.536 


0.0021 


2.8 


0.0025 


1.0 


0.0022 


3.6 


f'.f 6.2416 


0.349 


0.0027 


3.6 


0.0021 


3.6 


0.0020 


3.7 


ft 8.4184 


0.259 


0.0018 


3.0 


0.0016 


2.9 


0.0017 


3.8 


/(+ 3.8880 


0.560 


0.0022 


2.5 


0.0021 


2.9 


0.0018 


2.0 


f'+ 6.0648 


0.359 


0.0018 


2.2 


0.0022 


3.3 


0.0023 


3.3 


fg + 8.2416 


0.264 


0.0023 


3.0 


0.0019 


3.1 


0.0016 


2.7 


f' 4 + 10.4184 


0.209 


0.0017 


2.3 


0.0017 


2.9 


0.0013 


2.5 



3.2.1 The additional frequency of XY And 

We accept the second of the series listed in Table [3] 
(/2 = 7.7202) as an independent, additional frequency, 
because it has the highest amplitude both in the V and 
Ic bands. This frequency may be identified as a non- 
radial mode, as model calculatio ns (|Buchler fc Koilathll200ll ; 
iKollath. Molnar fc Sza bl l201ll) indicate that neither nor- 
mal nor 'strange'-modes have positive growth-rates with 
frequency ratio around 0.325, which is between the regime 
of the fifth and the sixth radia l overtones (see figure 3. in 
IKollath. Molnar fc Szabcfalll ). 

The explanation that these peaks belong to a secondary 
modulation series with ~ 0.196 cd -1 separation from the 
harmonics of the pulsation cannot be excluded either. How- 
ever, both the corresponding extremely short modulation 
period (~ 5.01 d) and the lack of the triplet components on 
the other side of the pulsation frequency make this explana- 
tion less likely. 



3.2.2 The additional frequency of UZ Vir 

In the case of UZ Vir, any of the /(, f^ and frequencies 
can be considered as the base-frequency of the additional 
frequency series; f[ and f' 3 are strong in the B band, while 
the amplitude of f 2 is the highest in the V and Ic bands. 
When fitting the final frequency solution to the light curves, 
fz was treated as a mathematically independent frequency, 
but this technical choice does not affect our subsequent re- 
sults and conclusions. 

As Fig. [3] shows, the ±lcd~ alias frequencies of the 
additional-frequency series of UZ Vir have similar ampli- 
tudes as the selected ones, which poses further ambiguities. 
Fitting the light curve of UZ Vir with all the detected fre- 
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quencies and with the series of listed in the first block of 
Table [3] the rms values of the residuals are 0.0112, 0.0093 
and 0.0096 mag in the 73, V and 7c bands, respectively. The 
solutions including the f( + and f[~ alias frequency series 
yield marginally higher, 0.0113, 0.0095, 0.0097 and 0.0114, 
0.0094, 0.0097 mag rms of the residuals. Moreover, all the 
amplitudes of the fijfi, fs and f' A terms in the B, V and 7c 
bands are higher than the corresponding amplitudes of the 
possible alias solutions. These results support our choice of 
the true additional-frequency series of UZ Vir. 

The frequency at 2.888 cd -1 (/{) has a period ratio of 
0.754 with fa. Although this is significantly larger than the 
observed fo j f\ period ratios of double-mode RR Lyrae stars, 
an open question is, whether in special, resonant cases, a 
low-amplitude, first-overtone mode might really be excited 
at such a large frequency ratio (Z. Kollath, private comm.). 
The frequency of the other possible base-component of this 
additional frequency series (j2 = 5.065 or = 7.242) might 
be identified as the third/fourth and the sixth-order radial 
mode or a non-radial mode. However, as RR Lyr models 
show that radial modes between the third and eights ra- 
dial orders are strongly damped, therefore this explanation 
seems to be unlikely. It is i mportant to note here al s o that , 
according to the study of iDziembowski fc Cassisil l)l999h , 
non-radial modes have the highest chances for excitation 
at close proximity of the radial modes. The explanation of 
the additional-frequency series as a secondary modulation is 
rejected in this case, because the peaks are at a distance of 
0.711 cd _1 from the pulsation harmonics, which corresponds 
to an unlikely short modulation period of only ~ 1.41 d. 

3.3 The final frequency solution 

None of the residual spectra shows any further significant 
peak after the removal of the frequency components iden- 
tified in the previous subsections. No frequency component 
with amplitude larger than 3.5u in at least one band is de- 
tected within the 0.05 d~ vicinity of the positions of the 
half integer pulsation frequencies (fc/o/2, k = 1,3, 5, 7, 9) in 
any of the two stars. 

The final frequency solution of XY And includes 57 
components: 13 harmonics of the pulsation frequency (fc/o), 
the modulation frequency (/ m ), 24 linear-combination fre- 
quencies belonging to the triplets (fc/o ± f m ), 15 quintuplet 
components (fc/o±2/ m ), the additional frequency (f? in Ta- 
ble O and 3 linear combinations of the additional frequency 
and the pulsation frequency. The rms values of the residuals 
are 0.011 and 0.008 mag in the V and 7c bands, respectively. 

As for UZ Vir, the total number of the identified fre- 
quency components is 84, of which 25 are harmonics of the 
pulsation frequency, one is the modulation frequency, 42 are 
triplet components, 11 are quintuplet components, one is a 
septuplet peak (6/o — 3/ m ), another one is the additional fre- 
quency, and three linear combinations of the additional fre- 
quency and the pulsation frequency are also included. The 
rms values of the residuals are 0.011, 0.009 and 0.010 mag 
in the B, V and 7c bands, respectively. 

The accepted values of the independent frequencies (fo, 
fm, fs) are determined by a non-linear Fourier fit to the 
V light curves of the stars, including all the identified lin- 
ear combination components with locked frequencies. These 
frequency values, the corresponding periods and their errors 



Table 4. Pulsation, modulation and additional frequencies, the 
corresponding periods and their errors (in parentheses, in the unit 
of the last digit) of XY And and UZ Vir. 



Frequency (cd 1 ) Period (d) 
XY And 


pulsation (fo, Po) 
modulation (/ m , 7m) 
additional (f 2 , P^) 


2.5079961(6) 

0.02417(3) 

7.7202(6) 


0.3987247(1) 
41.37(5) 
0.12953(1) 


UZ Vir 


pulsation (fo, Po) 
modulation (/ m , 7m) 
additional (f 2 , P^) 


2.1767879(4) 

0.014654(3) 

5.06481(6) 


0.4593925(1) 
68.24(1) 
0.197441(2) 



are listed in Table [4] The tables containing the full light- 
curve solutions (the frequency components, amplitudes and 
phases in each band) are available as supporting informa- 
tion with the electronic version of this article. Tables [5] and 
[6] are samples of these electronic data and give information 
regarding their form and content. 

The times of the pulsation and Blazhko maxima, that 
is, the zero pulsation and modulation phases, are determined 
from the light-curve solutions. The epoch of the pulsation 
corresponds to the time of the maximum of the mean pul- 
sation light curve, which is constructed as a synthetic light 
curve using the Fourier components of the pulsation fre- 
quency and its harmonics of the full light-curve solution. 
The epoch of the modulation is defined as the time when 
the amplitude of the fundamental frequency is the highest. 
This epoch does not necessarily coincide with the time of the 
highest amplitude phase of the modulation. The elements 
are given in the following equations: 



HJD 2 454 381.5513 + EpJ ls And ■ 0.3987247 d (1) 



Tgi mti ld = HJD 2 454 379.9 + E m ^ d And ■ 41.37 d 
TmlZ" = HJD 2 454 470.2294 + E^iT" • 0.4593925 d 



C = HJD 2 454 472.9 + E m z J iT 



68.24 d 



C2) 
(3) 
(4) 



4 COMPARISON AND DETAILS OF THE 
MODULATION OF XY AND AND UZ VIR 

As Fig. [T]shows, both the phase and the amplitude modula- 
tions of XY And are pronounced, while the fix point on the 
rising branch of the light curve of UZ Vir indicates that am- 
plitude modulation is dominant in this star. The phase of the 
maximum brightness of both stars varies, however, within 
about the same, ±0.02 pulsation-phase range in both stars, 
as can be seen in Fig. [4] This contradictory result warns that 
the definition of the amplitude of the phase modulation is 
not at all straightforward. 

The modulation of UZ Vir is of a rare type, because 
the amplitude of the low-frequency modulation components 
of the triplets are much larger than the high-frequency 
ones. The amplitude difference between the two modula- 
tion components of a triplet is characterized by the asym- 
metry parameter: Q = (A+—A-)/(A++ A-), introduced by 
lAlcock et all (|2003h . Both the MACHO (lAlcock et al.| | 2003 ) 
and the Konkoly Blazhko Survey I and II ( Sodor et alj|2012 ) 
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Table 5. Frequencies, amplitudes and phases of the V and Iq light-curve solutions (epoch = 54390.0, sine decomposition) of XY And. 
The full table is available as supporting information with the online version of this paper. 



XY And 








V 




Io 




identification 


/ 


amplitude 


phase 


amplitude 


phase 






fed - !> 


(mag) 


(rad) 


(mag) 


(rad) 


fn 




2.507996 


0.4263(3) 


5.275(1) 


0.2574(3) 


5.139(1) 


2 fn 




5.015992 


0.2044(3) 


0.275(2) 


0.1267(3) 


0.241(2) 


fm 




0.024171 


0.0101(3) 


3.11(3) 


0.0062(3) 


3.34(4) 


fo + 


fm 


2.532167 


0.0528(3) 


6.107(6) 


0.0320(3) 


6.127(8) 


2/o- 


1- fm 


5.040163 


0.0506(3) 


1.190(7) 


0.0317(3) 


1.278(8) 


fo- 


fm 


2.483825 


0.0225(3) 


2.07(2) 


0.0138(3) 


2.19(2) 


2/o- 


- fm 


4.991821 


0.0239(3) 


3.47(1) 


0.0150(3) 


3.53(2) 


3/o- 


h 2/ m 


7.572330 


0.0018(3) 


2.9(2) 


0.0016(3) 


3.0() 


4/o- 


t-2/ m 


10.08033 


0.0016(3) 


4.9(2) 


0.0009(3) 


5.6(3) 


fo- 


2/m 


2.459654 


0.0030(3) 


5.5(1) 


0.0019(3) 


5.8(1) 


2/o- 


-2/m 


4.967650 


0.0015(4) 


5.9(2) 


0.0010(3) 


5.2(3) 


/( = 


f'l ~ fo 


5.212223 


0.0026(4) 


5.5(1) 


0.0012(3) 


0.1(2) 


/ 2 




7.720219 


0.0029(4) 


0.4(1) 


0.0019(3) 


0.8(2) 


f* = 


fa + fo 


10.22822 


0.0021(4) 


2.7(2) 


0.0006(3) 


2.2(5) 


fi = 


/a + 2/o 


12.73621 


0.0017(4) 


4.0(2) 


0.0009(3) 


4.2(3) 



Table 6. Frequencies, amplitudes and phases of the B, V and Iq light-curve solutions (epoch = 54404.0, sine decomposition) of UZ Vir. 
The full table is available as supporting information with the online version of this paper. 



UZ Vir 









B 




V 




Ic 




identification 


/ 


amplitude 


phase 


amplitude 


phase 


amplitude 


phase 






(cd- 1 ) 


(mag) 


(rad) 


(mag) 


(rad) 


(mag) 


(rad) 


fo 




2.176788 


0.5485(4) 


0.867(1) 


0.4193(4) 


0.826(1) 


0.2603(4) 


0.681(1) 
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Figure 5. The amplitude vs phase variations (with uncertainties) of the first nine harmonic orders of the pulsation light curves of 
XY And and UZ Vir during their modulation cycles. Data of the bins nearest to 0.0 Blazhko phase arc marked with full dots. The solid 
lines show the same loops according to synthetic data. The amplitudes and phases are scaled by l k', the harmonic order, for the better 
visibility. Note that the phase range is the same in each panel: 0.6 rad. The horizontal axes are reversed in order to be in accordance 
with the variation of the phase of the maxima shown in Fig. [4] The directions of going around the curves are marked with arrows in the 
top-left panels. This direction remains the same for all the harmonic orders in both stars. 
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Figure 4. Brightness and phase variations of the pulsation-light 
maxima of XY And and UZ Vir in V band. Large (red) circles 
mark the observations in 11 and 10 Blazhko phase bins. Filled 
circles represent the bins nearest to zero Blazhko phase accord- 
ing to the elements given in Eqs. [2] and [4] Small (blue) dots are 
generated from synthetic light-curve solutions using the frequen- 
cies given in Tables [5] and \E\ Relative phases and magnitudes of 
the maxima with respect to the maximum of the mean pulsation 
light curve (marked by plus signs) are plotted. Arrows mark the 
directions of going around the loops; this is anti-clockwise for XY 
And and clockwise for UZ Vir. 



results show that the majority of the Blazhko stars have a 
positive asymmetry parameter. The asymmetry parameter 
of UZ Vir, calculated for the triplet around fo in V band, 
has an extremely large negative value: Q — —0.53 ± 0.01. 
Among the 731 Blazhko RRab stars in the MACHO LMC 
sample, only about 4 per cent ha ve larger negative Q values 
(see fig. 10 of lAlcock et~aLll2003h than this value. 

As ISzeidl fc Jurcsikl (|2009l ) pointed out, a connection 
exists between the sign of the Q parameter and the trav- 
eling direction of the light maximum along the loop in the 
maximum phase - brightness plane. The negative Q param- 
eter of UZ Vir reflects the property that its light maximum 
goes along the loop in clockwise direction (see Fig. 3]), while 
the maximum of most of the Blazhko stars either goes in the 
opposite, counter-clockwise direction (as XY And does), or 
the loop is degenerate. 

Examining Fig.[T] another significant difference between 
the modulation of the two Blazhko stars is discernible. The 
pulsation light curve of XY And is much smoother in each 
Blazhko phase than that of UZ Vir, especially around light 
minimum. The bumps and humps, characterizing the light 
variation of UZ Vir, are completely missing from the light 
curve of XY And. This difference is also reflected in the 
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Figure 6. Amplitude and phase relations of the triplet components of XY And and UZ Vir in the different harmonic orders. Squares 
(blue), circles (green) and triangles (red) denote the results in the B, V and Iq bands, respectively. Results for the different photometric 
bands of the same harmonic order are shifted horizontally for a better visibility. From top to bottom: (a) R^ = A kf +f m / A hfo (empty) 
and R^ = f Q _ f m / A^ f (filled) amplitude ratios of the triplet components relative to the harmonic components of the pulsation; 
(b) Q k = { A kf +f m - A klo-f m )/( A kf +fm + A fc/o-/m)> thc asymmetry parameter; (c) \t = <Pkf +f m ~ ¥>fc/ ~ <Pf m (empty) and 
X^T = <Pkf —f m ~ Vfe/o + ffm ( nuc( i), the phase differences between the modulation components and thc harmonics of the pulsation; 
(d) <£ fc = Vkf +f m + <Pkfa-f m - 2l Pkf Q (filled) and t> k = <Pkf a +f m ~ <Pkf -f m ~ 2< ?fm C^Pty) phase relations of the pulsation and the 
modulation components. 



light-curve solutions; although the noise characteristics of 
the residual spectra of the two objects are similar, signif- 
icantly more frequency components are needed (84) to de- 
scribe the more complex brightness variation of UZ Vir than 
for XY And (57). The harmonic components of the pulsa- 
tion are detected up to 25/o in UZ Vir, while only up to 
13/o in XY And, which indicates that the mean pulsation 
light curve of XY And is more simple than that of UZ Vir, 
as well. 



4.1 Amplitude and phase relations of the 
pulsation light curve 

Fig. [5] shows the amplitude vs phase variations of the first 
nine harmonic orders of the pulsation light curves of XY And 



and UZ Vir in 11 and 10 phase bins of their modulation 
cycles, respectively. The solid lines show the same loops de- 
rived from synthetic light-curve solutions. The directions of 
going around the curves are marked with arrows in the top- 
left panels. This direction remains the same for all the har- 
monic orders in both stars. 

The phase - amplitude loops are dissimilar in the differ- 
ent harmonic orders for both stars and none of the loops 
matches exactly the observed maximum phase - maximum 
brightness loops shown in Fig. 2] Furthermore, the max- 
imum amplitude of the loops occurs in different Blazhko 
phases, as can be seen in Fig. [5] where the phases and am- 
plitudes of the light curve at zero Blazhko phase are denoted 
by solid symbols. Some of the highest order loops show an 
interesting feature; they seem to 'open' at the lowest am- 
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plitudes in the 9th, and the 8 - 9th orders for XY And and 
UZ Vir, respectively. In accordance with its more complex 
light variation, the phase - amplitude loops of UZ Vir show 
a large variety, from a circular shape (1st order) to a quasi- 
degenerate one (7th order), while, with the exception of the 
9th order, all the phase - amplitude loops of XY And are 
'egg'-shaped or elliptical. 

4.2 Amplitude and phase relations of the triplet 
components 

The relative amplitudes of the pulsation and modulation 
components, 

Rk = Akf +f m /Akf 0: (5) 
Rk = A kf0 - f JA kf0 , (6) 
and the asymmetry parameter 

Qk = (A kfo+fm - A kf0 _ fm )/{A kfo+fm + A fcJW J (7) 

are derived for each harmonic order (fc). 

To follow the phase variations, the epoch-independent 
phase differences between the pulsation and the modulation 



components, 

xt = fkfo+Jm - fkfo - <Pfm (8) 
Xk = <Pk} a -f m - VkJo + <Pfm ( 9 ) 

and phase relations of the two modulation components, 
*fe = <Pkf +f m +^fe/ -/ m - 2 ^fc/ (10) 
*fc = Vkfo+U - Vkf -f m - 2^/ m (11) 



are calculated according to the Fourier solutions of the light 
curves given in Tables [S] and [U 

Fig. [B] shows the R+, , Q k , xt , X k > *fc> and *fe am- 
plitude and phase relations of the triplet components in the 
different harmonic orders for the observed bands of the two 
studied stars. 

4.2.1 XY And 

Both the amplitude and phase variations show smooth, 
gradual, almost monotonic changes. The first panel (a) il- 
lustrates that while the low-frequency triplet components 
hardly changes as compared to the pulsation components, 
the amplitudes of the high-frequency components become as 
high as the amplitudes of the pulsation frequencies from the 
7th harmonic order. As a result, the asymmetry parameter, 
Q k , is increasing towards higher harmonic orders (panel b). 

The bottom set of panels displays the phase-difference 
variations of the triplet components. Each of the four phase 
relations shows quite a regular behaviour. It is interesting 
to note that while only 0.5 rad gradual increase of the phase 
difference between the higher amplitude triplets and the pul- 
sation components (xt) is observed, the phase difference of 
the lower amplitude triplet component (x k ) is less stable, 
it changes within about 1.5 rad. The $ fc and ty h phase rela- 
tions also show monotonic changes with increasing gradients 
towards higher harmonic orders. 

The most interesting feature of the phase relations is 
that each of the four shows a definite colour dependence, 
which is the most pronounced for x k and ^fe- These phase 



relations in the V and 7c bands differ systematically in each 
harmonic order by 0.3 -1.0 rad. The Q k asymmetry parame- 
ter may also be somewhat colour dependent; the asymmetry 
of the triplets seems to be slightly larger in the Ic band than 
in the V band. This tendency, however, is not really signifi- 
cant. 



4.2.2 UZ Vir 

The amplitude and phase relations of the triplets of UZ Vir 
(right-hand panels of Fig. [6} show a more complex behaviour 
than for XY And. Most probably, this is the consequence 
of the appearance and disappearance of the bump/hump 
features on the light curve during the modulation cycle. 

The amplitude variation of the lower amplitude triplet 
component, Rt , is smoother than that of the higher ampli- 
tude triplet component, , which shows local maxima at 
the 5th and the 9 - 10th orders. The variation of the asym- 
metry parameter is also complex, it has local minima (tak- 
ing into account the negative values of Q k , these correspond 
to large asymmetry) in the first, the 4 -5th and 9 -10th 
harmonic orders. Towards higher orders, the asymmetry of 
the triplets (\Q k \) shows a decreasing tendency, while it in- 
creases for XY And. 

The phase behaviour of the triplet components of 
UZ Vir is not smooth either. The total ranges of the varia- 
tions of xt i Xfe an d & k are 1.5-2 rad, while the variation 
of *& fe extends to 2.5 rad. The x k phase difference follows 
similar changes to the variations of the asymmetry param- 
eter with local minima in the 3rd and 8th orders. With the 
exception of the first order of $ fe , the variations of and 
$ fe reflect the changes of xt an d x k > respectively. 

Colour dependence of the phase relations of UZ Vir 
is also detected. The xt an d ^ k values are systematically 
larger by 0.1-0.5 rad in the Ic band than in the B and 
V bands, while the colour dependence seen in the x k data 
is the opposite. With the exception of the first order, the 
asymmetry parameter of UZ Vir is also colour dependent. 
Contrary to the results for XY And, the triplets are less 
asymmetric towards longer wavelengths here. 



5 PHYSICAL PARAMETERS AND THEIR 
VARIATION DURING THE MODULATION 

The variations of the pulsation-averaged atmospheric pa- 
rameters of XY And and UZ Vir during the Blazhko cycle 
are derived us ing the inverse photometric B aade-Wesselink 
method (IPM; ISodor. Jurcsik fc Szeidj|200gh . These param- 
eters are the effective temperature (T e g), luminosity (L), ra- 
dius (R) and effective surface gravity (logg). Utilizing mul- 
ticolour light cu rves and synthetic colour s from static atmo- 
sphere models (|Castelli fc Kuru"c3l2003l ). the method finds 
the best fitting model parameters through the variation of 
template T a a and U ra d curves. In order to test the stability 
and reliability of the method, different V rll d template curves 
and different weightings are applied. Therefore, similarly to 
the previous analyses, the IPM cod e is run with four differ- 
ent in ternal settings (for details see lSodor. Jurcsik fc Szeidll 
120091 . Table 1) to estimate the inherent, method-specific un- 
certainty of the results. 
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Table 7. Mean physical parameters and their uncertainties derived from the mean pulsation light curves of XY And and UZ Vir with the 
inverse photometric method (IPM) using static atmosphere models with different [Fe/H] metallicities. The luminosities, temperatures 
and masses are constrained to match the evolutionary possible values for horizontal-branch stars. The mean physical parameters are also 
required to fit the pulsation period via the pulsation equation. (My) denotes the intensity-averaged absolute V magnitude. 



Object [Fe/H] W d (M v ) R L T cB 

(dex) (OT ) (pc) (mag) (R ) (L ) (K) 



XY And 


-0.68 


0.58 ±0.01 


3700 ± 200 


0.76 ± 0.05 


4.15 ±0.20 


40.5 ±2.5 


7100 ± 100 


UZ Vir 


-0.65 


0.61 ±0.01 


2850 ± 150 


0.77 ±0.05 


4.54 ±0.20 


41.0 ±2.5 


6800 ± 100 


UZ Vir 


-0.90 


0.59 ±0.01 


2900 ± 150 


0.71 ± 0.05 


4.55 ±0.20 


43.5 ± 2.5 


6900 ± 100 


UZ Vir 


-1.26 


0.57 ±0.01 


2950 ± 150 


0.65 ± 0.05 


4.56 ±0.20 


46.0 ±2.5 


7000 ± 100 



5.1 Constant parameters 

To apply the IPM, the Blazhko-phase-independent param- 
eters have to be known first. These are the metallicity 
([Fe/H]), mass (9JZ), distance (d), and the mean dereddened 
colours of the variables averaged over both the pulsation and 
modulation periods. 

Th e metallicity o f XY And was me asured by spectro- 
scopi c (|Lavdenl Il994l ) and photometric (|Jurcsik fc Kovacsl 
1 19961) means. Both met hods yield the same, —0.68 dex re- 
sult (jjurcsik et al.|[2009l . table 5). Therefore, this value, as a 
reliable estimate for the metallicity of XY And, is accepted 
and applied throughout the following analysis. No spectro- 
scopic measurement of the metallicity of UZ V ir has been 
published. It was shown bv lJurcsik et al.l (|2009T l that a well 
defined mean pulsation light curve of a Blazhko star yields a 
reliabl e metallicity value accor ding to the [Fe/H](P, $31) for- 
mula (|Jurcsik fc Kov acs 1996) even for strongly modulated 
RR Lyrae stars. This method gives [Fe/H] = —0.90 dex for 
UZ Vir. However, the example of RZ Lyr warns that, in cer- 
tain cases, the photometric metallicity may di ffer by as large 
as ~ 0.4 dex from the spectroscopic value (|jurcsik et al.l 
I2012T I. Therefore, the effect of the choice of the input metal- 
licity on the IPM results is checked by calculating the mean 
physical parameters of UZ Vir using two other metallicity 
values, as well. As only evolutionary possible values of the 
mean physical parameters are accepted, the other metallicity 
values (—0.65 and —1.26) were selected from the [F e/H] val- 
ues of the horizontal branch evolutionary models of lDormanl 
l|l992T ). 

Our differential photometric observations do not allow 
the determination of the mean (reddened) colours of these 
two variables, since no standard multicolour magnitude mea- 
surement of the comparison st ars or of any other obj ect in 
the field of view is available. I Schmidt fc Sethi (|l996f ) give 
the intensity-averaged values of (V) — (R) = 0.22 ± 0.02 
and 0.24 ± 0.03 mag for XY And and UZ Vir, respectively. 
However, these colour indices should be treated with cau- 
tion, since they are calculated only from 10 - 20 data points 
per band, observed in different Blazhko phases. Therefore, 
the uncertainties of these colours are estimated to be a 
few hundredth magnit ude higher than the values given by 
Schmidt fc Sethi (Il996lh According to the extinction maps of 



Schlegel et all (|l998T l. the interstellar reddening [E(B - V)] 



towards the directions of XY And and UZ Vir are 0.050 and 
0.026 mag. Since these values also have uncertainties on the 
order of 0.01 mag, and because these are only upper limits 
for the actual reddenings, the possible colour ranges, con- 
sequently the possible overall temperature ranges, are too 



wide to impose any meaningful constraint on the IPM fit- 
ting process. 

Therefore, the temperature, luminosity and mass of 
both stars have been constrained using the pulsation equa- 
tion and horizontal-branch (HB) evolu tionary-model results . 



2003, 



1992, 



According to the pulsation equation of lMarcpni et al.l 
eq. la) and the HB evolutionary models of lDormanl 
figs. 7-10), the pulsation periods of the studied objects cor- 
respond to the following physical parameter ranges: 7000 - 
7200 K and 6700 - 7100 K temperature, 38 - 43 L@ and 38.5 - 
48.5 Lq luminosity for XY And and UZ Vir, respectively. 

Applying these constraints, the IPM is run for the mean 
pulsation light curves of both objects to determine consis- 
tent sets of the mean and constant parameters. The re- 
sults from these runs are listed in Table [7] To check the 
effect of the uncertainty in the metallicity of UZ Vir, the 
mean physical parameters are calculated for three possi- 
ble metallicity values. The di stances are estimated a ccepting 
(V) = 13.73 and 13.05 mag (|Schmidt fc Sethll 19961 ) and as- 
suming A(V) = 3.14E(B-V) = 0.157 and 0.082 interstellar 
absorption for XY And and UZ Vir, respectively. The un- 
certainties given in Table [7] are estimated by taking into 
account the inherent error of the method, the uncertain- 
ties of the input parameters and the possible ranges of the 
constrain ed parameters a ccording to the HB evolutionary 
models of lDormanl (|l992T ). The data given in Table [7] prove 
that the three possible metallicity values of UZ Vir results in 
similar uncertainties of the output parameters as the other 
factors together. The only exception is the radius, which is 
practically the same for all three input metallicities. 



5.2 Blazhko-phase-dependent parameters 

Fixing the Blazhko-phase-independent physical parameters 
to their values determined in Sect. l5.ll and listed in Table [TJ 
the IPM is run using the light curves of XY And and UZ Vir 
in 11 and 10 different, disjunct phase bins of the modulation, 
respectively. In order to obtain the highest possible number 
of Blazhko phase bins with full pulsation-phase coverage, 
the bins are not exactly equidistant, but spread over the 
Blazhko cycle almost evenly. Each bin is narrow enough for 
the pulsation light curves not to be affected by Blazhko- 
phase smearing. 

The IPM is run for UZ Vir with the three possi- 
ble metallicities listed in Table [7J to check its effect on 
the results. It is found that the different metallicity val- 
ues affect only the means of the physical parameters, but 
the amplitudes and the phase relations of their variations 
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Figure 7. The mean observed (left-hand panels) and the IPM-derivcd physical parameters (right-hand panels) of XY And and UZ Vir in 
different phases of the modulation. From top to bottom, the left-hand panels show: amplitude modulation - amplitude of the /o frequency 
component; pulsation-phase or period modulation - variation of the phase of the /o pulsation component (empty circles) and deviation 
of the instantaneous pulsation period from the mean pulsation period (continuous line); pulsation-averaged V magnitude; pulsation- 
averaged B — V colour; pulsation-averaged V — Iq colour. Magnitude- and intensity-averaged brightnesses and colours are denoted by 
angle and round brackets, respectively. From top to bottom, the right-hand panels show the pulsation averages of the following physical 
parameters: radius, surface gravity, absolute visual magnitude, luminosity, and effective temperature. Results of four different settings of 
the IPM are shown in order to indicate the inherent uncertainity of the method. 



with Blazhko phase remains practically the same. There- 
fore, only the results obtained with the most likely value of 
[Fe/H] = — 0.90dex are shown. 

The IPM gives the radius, luminosity, temperature and 
surface-gravity changes during the pulsation in each phase 
bin of the modulation. The pulsation-averaged mean values 
of these quantities characterize the global mean changes in 
the stellar parameters during the modulation cycle. Simi- 
larly to the IPM results of the previous investigations of 
Blazhko stars, the actual values of the fixed or constrained 
physical parameters influence the averages of the other phys- 
ical parameters, nevertheless, they affect the variations dur- 
ing the Blazhko cycle only slightly. The amplitudes of the 
solutions differ somewhat if the input parameters are varied 
within their possible ranges (Table [7J, but the relative sign 
of their variations remains unchanged. 

The appearance of f m in the frequency solution of the 
light curves reflects the variation of magnitude averaged 
mean brightness in the different bands during the modu- 
lation. The IP method, utilizing colour information as well, 
discloses what luminosity, radius and temperature changes 
are required to describe the observed changes in the mean 
magnitudes. The IPM results on XY And and UZ Vir are 
plotted against the modulation phase in Fig. Here, the ob- 
served (left-hand panels) and IPM-derived (right-hand pan- 
els) variations of the pulsation-averaged parameters during 



the Blazhko modulation are plotted. The observed param- 
eters (from top to bottom) are the amplitude of the main 
pulsation frequency (/o); the $1 phase variation of this fre- 
quency component and the pulsation-period variations de- 
duced from the observed phase variations of the /o compo- 
nent; and the different averages of the light and colour curves 
according to magnitude- and intensity-scale representations. 



6 DISCUSSION AND SUMMARY 
6.1 Light-curve modulation 

The Fourier spectra of the light curves of XY And and 
UZ Vir are described by multiplet frequencies (kfo ± n/ m ) 
and the modulation frequency (/ m ), typical of Blazhko stars. 
The multiplet structures are triplets and quintuplets, with 
a single septuplet component in the spectrum of UZ Vir 
(6/o — 3/ m ). The phase of the rising branch of XY And is 
strongly modulated, while a fix point on the rising branch 
of UZ Vir (see Fig. [TJ indicates that amplitude modulation 
dominates its light-cure variation. Looking at the phase vari- 
ations of the maximum brightness of the two stars, similar- 
amplitude changes are observed as shown in Fig. [3] This 
controversial behaviour of the phase variations of the rising 
branches and the maxima warns that the strength of the 
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phase modulation at different parts of the pulsation light 
curve can be rather different. 

The pulsation light curve of XY And is much smoother 
in each Blazhko phase than that of UZ Vir. Consequently, 
the light-curve solution of XY And involves less pulsation 
and modulation components than that of UZ Vir. 

The modulation of UZ Vir is of a rare type, since the 
amplitude of the low-frequency triplet components are sig- 
nificantly larger than that of the high-frequency ones. Ac- 
cordingly, the light maximum goes in clockwise direction 
around the loop on the pulsation phase - brightness plane 
(Fig. while the maximum of most of the Blazhko stars 
goes around in the opposite direction, or the loop is degen- 
erate. 

For the first time, the amplitude and phase relations of 
the triplet components are also studied in detail. No com- 
mon feature of the variations of the defined parameters of 
the two stars has been recognized. This result supports that 
Blazhko stars behave very individually. It is difficult, if at all 
possible, to find any overall, common property of the mod- 
ulation. The most interesting finding of this investigation 
is the detection of the colour dependency of the asymme- 
try parameter and the phase relations. The colour depen- 
dence of these quantities is the opposite of each other in the 
two stars, which might be connected to their opposite-sign 
amplitude and phase variations during the Blazhko cycle. 
The colour dependence of the amplitude and phase rela- 
tions of the triplets of other Blazhko stars has to be studied 
to decide whether such a connection really holds. As the 
detected colour dependence carries important physical in- 
formation on the modulation, its further investigation and 
explanation should be an important step in the description 
and interpretation of the Blazhko effect. 

6.2 Additional frequencies 

The additional frequencies in RR Lyrae stars have been dis- 
covered only recently, as very precise and extended photom- 
etry is needed to find these low-amplitude signals. 

Both XY And and UZ Vir show variability besides the 
pulsation and modulation with an additional, independent 
frequency. Series of four peaks appear in the Fourier spectra 
of both objects with separation that equals to their pulsation 
frequency. One of these peaks is identified as a linearly inde- 
pendent frequency component, while the other three mem- 
bers of the series are interpreted as linear combinations of 
this and the pulsation frequency. The frequency ratios of the 
possible main components of the additional frequency series 
(fo/f) are 0.325 for XY And and 0.754, 0.430 or 0.301 for 
UZ Vir. The 0.301, 0.325 and 0.430 frequency ratio values 
cover the possible regime of the third-sixth radial overtones. 
Unfortunately, the actual frequencies of the higher order ra- 
dial modes on a large grid of stellar parameters have never 
been published. Thus we do not know whether or not the 
observed frequency ratios fit any higher order radial mode, 
indeed. However, as RR Lyrae models indicate that these 
higher order modes are strongly damped, therefore we con- 
sider it unlikely that radial modes were detected in these 
stars. 

In the case of XY And, two possible explanations of 
the additional-frequency series have been raised: (a) one fre- 
quency of the series might be a non-radial mode, (b) all 



four frequencies correspond to a secondary modulation of 
an extremely short, ~ 5.01 d period. For explaining the se- 
ries of additional frequencies in UZ Vir, we also raised two 
possibilities: (a) the frequency at 2.89 cd _1 might be the 
first radial overtone, however, with an unusually large /o //i 
ratio of 0.754, (b) one of the two highest amplitude com- 
ponents of the additional frequency series corresponds to a 
non-radial mode. The recent detection of an additional fre- 
quency with 0.753-758 frequency ratio in the Kepler data of 
RR Lyr (Molnar et al. submitted to ApJ), which has been 
identified with the first overtone mode, favours the first pos- 
sibility in the case of UZ Vir, too. 

It is of interest to examine the reported occurrences 
of other, additional frequencies among RR Lyrae vari- 
ables. Both double mode RR Lyrae stars with extended 
satellite observations (AQ Leo and Co RoT ID 0101368812 
iGruberbauer et al.|[2007l : IChadidllioi^ ) show additional fre- 
quencies beyond the fundamental and the first overtone 
modes. Many Blazhko variables have one or more extra, in- 
dependent frequencies in their light variation in addition to 
the pulsation and modu lation compon e nts. S u ch frequencies 



have been reported by Chadid et al. ll201Cf) . [Poretti et al.l 



j2010t ) a nd ICuggenberger et alj (|201ll . [201 



ble 2. of iBenko et all l|2010D . The ±12.5/, 



see also Ta- 
eaks around 

some o f the pulsation components in MW Lyr (|jurcsik et alj 
l2008al ) might al so be connected to a n additional frequency, 
as suggested bv lPoretti et al.l ([2010T I. In contrast, only two 
of the 19 non-modulated RRab stars observed by Kepler 
show frequencies b esides the fundamental-mode pulsation 
(N emec et al]|201ll ). and no additional frequency has been 
detected in the only non-Blazh ko CoRoT RRab star anal- 
ysed so far ([Paparo et al.l 120091 ') . In the residual spectra of 
the 14 non-modulated RRab stars observed in the Kon koly 
Blazhko Survey I (see references in ljurcsik et al.ll2009l ). no 
additional frequency has been found either. Although the 
majority of the Kepler Blazhko targets still lack thorough 
analysis, it is clear that the occurrence rate of the appear- 
ance of additional frequencies is much higher among mod- 
ulated RRab stars than in non-modulated ones. A possible 
explanation might be that the changing physical parameters 
during the Blazhko-cycle ([Jurcsik et al.l l2008bT ) temporar- 
ily fulfill the conditions required for the excitation of the 
other observed (radial or non-radial) modes, whereas the 
chance to satisfy the excitation conditions is much smaller 
for stable-light-curve RRab stars. A similar explanation has 
been suggested for the appearance of a peculiar bump on the 
descending bran ch of RZ Lyr in the low-amplitude phases of 
its modulation (|jurcsik et al.l [20121 ). However, the excited, 
additional mode of RZ Lyr is supposed to be in resonance 
with one of the harmonics of the main pulsation frequency, 
which explains why it does not appear in the Fourier spec- 
trum as an additional frequency directly. 



6.3 Physical-parameter changes during the 
modulation 

The modulational changes in the atmospheric parameters 
of XY And and UZ Vir have been derived applying the 
IP method. Contrary to the many differences between the 
Blazhko-modulation properties of the two objects, their 
mean physical parameters are rather similar; XY And is only 
~150K hotter, ~4L© fainter, ~0.4_R Q smaller, ~0.01 OT Q 



The Blazhko modulation of XY And and UZ Vir 13 



less massive and ~0.22 dex more metal rich, than UZ Vir. 
The changes in the mean temperature, radius and luminos- 
ity during the Blazhko cycle of the two stars do not show 
any significant difference either. 

Both objects show strong luminosity variation during 
the modulation: AL « 2 Lq and » 2.5 Lq are derived for 
XY And and UZ Vir, respectively. Similarly to all Blazhko 
RR Lyrae stars analysed with the 1PM so far, XY And and 
UZ Vir are the brightest when the pulsation amplitude is 
the highest. 

The radius variation is weak; the average radius is the 
highest around Blazhko maximum in both variables. The 
IPM results show only marginal temperature variation of 
XY And during the Blazhko cycle, but the T e g and (V) — 
(Ic) curves together suggest that XY And is about 10 - 20 K 
warmer at the high-amplitude phase of the modulation than 
when the pulsation amplitude is the lowest. The IPM results 
on UZ Vir indicate only a modest but definite modulational 
variation of about 30 K in temperature. Both the radius and 
temperature are the highest at Blazhko maximum both in 
XY And and in UZ Vir. 

The Blazhko model proposed bv lStothers! (|2006l ) postu- 
lates a changing magnetic field in the envelope of the star, 
which influences the parameters of convection. The model 
makes predictions about the phase relation between the 
amplitude and period modulations and between t he period 
modulation and radius variation (|Stothers! l201ll ). namely, 
these relations dep end on the temp erature of the star. For 
hot RRab stars, the lStothersI (|2006h model predicts that the 
pulsation period is decreasing at high-amplitude Blazhko 
phase (the period and amplitude variations are in anti- 
phase), while for stars below the crossover temperature of 
about 6400 K, the period is increasing at the same Blazhko 
phase (the changes of the two quantities are in phase) . These 
relations are shown in the top two panels in the left columns 
of Fig. [7] for the two studied objects (the pulsation-period 
changes are plotted with continuous line). The two vari- 
ables have similar temperatures well above the predicted 
crossover value, and they also have similar pulsation periods 
and metallicities. Still, the phase relations between their pul- 
sation amplitude and pulsation period variation are the op- 
posite of each other. This behaviour is also demonstrated in 
Fig. [4] which shows that maxima of the two objects go along 
their modulational loops in the pulsation phase - brightness 
plane in opposite direction. 

The Stothers model also expects a relation between the 
pul sation period a nd radius variation during the Blazhko cy- 
cle (|Stotherdl20~ilT ). The model predicts anti-correlated vari- 
ation between these two quantities in hot RRab stars, which 
turns into correlation for cooler RRab variables. However, 
both objects are hot, short-period RRab stars, and both 
have the highest radius at Blazhko maximum, while their 
pulsation-period variations are the opposite of each other. 

These findings are not in accord wit h the predictions of 
the Stothers model (|Stothersll200f3 . l201ll '). which claims that 
two Blazhko stars with such similar properties are expected 
to display similar relations between their amplitude, pul- 
sation period and radius variations. Our results show that 
even if the fundamental idea of the Stothers model is right, 
the relations between the studied properties of the Blazhko 
stars are more complicated than predicted by the model in 
its present form. 
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Table Al. CCD AV observations and derived colour in- Table A4. CCD AB time series of UZ Vir. 

dices of XY And, relative to the comparison star 2MASS 

01270016+3404307. HJD - 2400000 AB(mag) 



HJD - 2400000 


AV(mag) 


AV — /(mag) 


54504.57158 
54504.57737 


0.975 
0.986 


54390.46709 


0.466 


-0.211 






54390.47001 


0.483 


-0.201 







Table A5. CCD A7 C time series of UZ Vir. 



Table A2. CCD AJ C time series of XY And. 



HJD - 2400000 A7c (mag) 



54390.46855 0.677 
54390.47148 0.678 



HJD - 2400000 A/ c (mag) 



54504.57515 1.309 
54504.58108 1.296 



APPENDIX A: SAMPLES OF THE 
ELECTRONIC TABLES 

The full tables, containing the observations of UZ Vir and 
XY And (Tables [A1HA5[1 . are available as supporting infor- 
mation. 



Table A3. CCD AV observations and derived colour in- 
dices of UZ Vir, relative to the comparison star 2MASS 
13082756+1322403. 



HJD -2400000 AV(mag) AB - V(mag) AV - /(mag) 



54504.57397 1.183 -0.210 -0.109 

54504.57976 1.191 -0.217 -0.102 



